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1. Introduction 



There are several experimental evidences on power-law tailed statistical distribu- 
tions of single particle energy. One can find power-law tailed distributions also in 
a wide range of other fields of science [ ffl Gl CU ■ Our aim here is to reproduce 
exponentially starting and power-law tailed spectra in the framework of the non- 
extensive thermodynamics, instead of fitting the low- and high-energy parts of the 
spectra separately. To achieve our goal we use a generalized Boltzmann equation [ 
U) which accomodates a general, associative energy composition rule, defining this 
way a certain realization of non-extensive thermodynamics. Another way to gener- 
alize the Boltzmann equation preceeded our approach, it can be connected to 
our extension of the Boltzmann theory by requiring the same stationary solution [ 
El- To each composition rule a strict monotonic function can be constructed which 
transforms it to a normal extensive quantity. In the theory of formal groups this 
function is called the logarithm of the group. We investigate also the evolution of 
the one-particle distribution in the framework of this two-body Boltzmann equation 
supported by a non-extensive energy composition rule. 
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2. Non-Extensive Boltzmann Equation 

In the non-extensive extension of the Boltzmann equation we keep the original 
factorizing form for the two-body distributions, 

d f 

ol/l = / W1234 (/a/4 - /1/2) , (1) 
OT J 234 

but in the transition rate 

^1234 = M? 234 ■ S ((ft + p 2 ) - (p 3 + &)) • <5 (/i (Ei,E 2 ) ~ h (E 3 ,E 4 )) , (2) 

the two-body energy composition is not necessarily extensive, {En = h (Ei,E 2 ) 7^ 
Ei + E 2 ). 

Clearly Ei — \Jp\ + irif are regarded as asymptotic, purely kinetic energies, 
while h(Ei,E 2 ) may contain further contributions stemming from the pair inter- 
action. It is not trivial whether these can always be divided to one-particle con- 
tributions, supporting a quasi-particle picture. In the followings we demonstrate 
that under quite general assumptions about the value of the function h(x,y) it can, 
however, be done. 

We assume that the generalized energy sum is associative 

h(h(x,y) ,z) = h(x,h(y,z)); (3) 

a reasonable assumption for approaching the thermodynamical limit. Due to a 
mathematical theorem one finds a function X(h), which maps the energy composi- 
tion rule to additivity, i. e. 

X(h)=X(x)+X(y) (4) 

is fulfilled. This solution is unique up to a constant multiplicative factor. In this 
case the stationary solution of the non-extensive Boltzmann equation is given by 

/(p) = Ie-^)/ T . (5) 

Table ^ shows a couple of interesting composition rules, which lead to widely 
studied statistics. 

3. Numerical simulation 

We solve numerically the non-extensive Boltzmann equation by a parton cascade 
simulation. We use several different initial momentum distributions, and then make 
random binary collisions between randomly chosen pairs of particles. By doing so 
we apply the following rules: 



X(Ek)+X (E 2 ) = X (E 3 ) + X (£4) 



(G) 
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h(x,y) 


X(E) 


Z ■ f(E) 


Gibbs 


x + y 


E 


e -E/T 


Tsallis 


x + y + axy 


Mn(l + aE) 


(1 + aE)- 1 /^ 


Levy 


(x b + y b )^ b 




e ~(aE) b /{aT) 


Renyi 


axy 




{aE)- 1 /^ 



Table 1. Examples for the h(x,y) composition rule. 



Pi + P2 = P3 + Pi (7) 

After selecting two colliding particles in the simulation, we find the value for the new 
momentum of the first particle (pz) satisfying the above constraints but otherwise 
random. Then applying the rules in eqs. © and J7J) we can calculate the momentum 
of the second outgoing particle {pa)- Presently we use the free dispersion relation for 
massless particles (Ei(pi) — \pi\). A typical simulation includes 10 6 — 10 7 collisions. 
After 3—5 collisions per particle on the average, the one-particle distribution reaches 
its stationary form. 

There are several conserved quantities in the simulation, 

N N 

X (E tot ) = J2 X ( E i) . p = X>> N 

i=l i=l 

which are monitored during the simulation. 

4. Kinematics with non-extensive rules 

Most power-law tailed spectra emerging from the Tsallis thermodynamics [ [7] are 
based on the following composition rule 

h(x, y) = x + y + axy. (9) 

It can be regarded as the sub-leading order expansion of eq. f° r small arguments. 
The traditional Boltzmann distribution is recovered in the case of a = 0. 

In one certain binary collision the total quasi-energy of the two particles is 
conserved, so after choosing randomly a new energy value for the first outgoing 
particle (E3), we can calculate the new energy value of the second outgoing particle 
(E4), by inverting eq. 10, 

£4= E ™~!f 3 , E 12 =E 1 +E 2 + aE 1 E 2 . (10) 
1 + aiLi, 



N 

"E 1 ' (8) 

i=i 
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Fig. 1. Kinematics with h(x. y) = x + y + axy rule for E = 4 and P = 3.5. 



On Fig. ^a- 11 example is shown, how one can obtain p^ and p*3 when the energy 
and the momenta of the incoming two particles are known. This figure shows cross 
sections of the surfaces in the phase space, whereon the endpoint of p\, which is 
also the startpoint of P2, lies on. p% +p^ = p\ +P2 = P, and the new momentum p% 
should also point to this surface. For the composition rule in eq. © the endpoints 
of the vector p^ constitute a surface of revolution of fourth-order curves. The special 
case of an ellipse emerges for a = 0. 

5. Spectra with non-extensive rules 

The result of three simulations with the energy composition rule of eq. © is shown 
on Fig. [3 In all of the three cases the number of particles (10 5 ) and the number 
of collisions (10 6 ) are the same, the particles are massless. Eg is chosen to get the 
same value for X(E to t) for each value of the parameter a. 

The a = case leads to the traditional Boltzmann distribution (curve 2) , while 
the one with positive Tsallis parameter (a = 0.5) has a power-law tail above the 
exponential curve (curve 3) and the distribution with negative parameter (a = 
—0.25) goes below the exponential (curve 1). 

6. Conclusions 

We demonstrated how to simulate power-law tailed stationary distributions in a 
parton cascade code. As a generic example the cut power-law distribution was 
generated. We have found a unique non-exponential stationary distribution for 
each non-extensive energy addition rule. The non-extensivity parameter is yet to 
be derived from microscopic interactions. We conclude that power-law tails in one- 
particle spectra are compatible with stationary states of generalized non-extensive 
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Energy distributions with Tsallis rules 

1. : a = -0.25, E = 0.294 GeV 

2. : a = 0.00, E = 0.350 GeV 

3. : a = 0.50, E = 0.490 GeV 



■Tr^2x(Ei) = 0.987 GeV 



0.01 0.1 1 10 

E [GeV] 

Fig. 2. One-particle energy distributions with h(x, y) = x + y + axy rule. 

kinetic models. 
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